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Abstract. Future astrometric satellites, such as SIM (NASA's Space 
Interferometric Mission) and GAIA (ESA's Global Astrometric Interfer- 
ometer for Astrophysics), hold the promise of mapping out the detailed 
phase space structure of the Galactic halo by providing unprecedented 
annual proper motion and parallax of 1 — 10/ias astrometric accuracy. 
Here we show that proper motions of hundred or so giant branch stars 
in a tidal debris torn from a small satellite (a 10 5 Lq Galactic dwarf 
galaxy or globular cluster) in the halo is sensitive to the current Galactic 
potential and its past evolution. We follow the evolution of a cold (veloc- 
ity dispersion of 10 km/s) stream on a nearby (between 8-50 kpc) polar 
orbit in a variety of histories of the potential of the Galaxy, and observe 
the bright (V < 18mag) members of the debris tail with GAIA accuracy. 
We simulate effects due to the growing or flipping of the Galactic disk 
over the past 4 Gyrs or the perturbation from a massive accreted lump 
such as the progenitor of the Magellanic Clouds. We study various factors 
influencing our ability to identify streams, including contamination from 
field stars, accuracy of radial velocity and distance data and evolution 
and non-axial symmetry of the potential. Our simulations suggest that 
nearby, cold streams could be detected with GAIA if these cousins of the 
Sagittarius stream exist. Results of Johnston, Zhao, Spergel & Hernquist 
(1999) and Helmi, Zhao & de Zeeuw (1999) for static Galactic potentials 
are likely to be largely generalizable to moderately time-evolving poten- 
tials. SIM and GAIA measurements of debris stars might be used to 
probe both Galactic structure and Galactic history. 



1. Introduction 

Tidal streams in the Galactic halo are a natural prediction of hierarchical galaxy 
formation, where the Galaxy builds up its mass by accreting smaller satellite 
galaxies. They are often traced by luminous horizontal and giant branch stars 
outside the tidal radius of a satellite (such as the Sagittarius dwarf galaxy) or 
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a globular cluster (cf. Grillmair et al. 1998, Irwin & Hatzidimitriou 1995) as 
a result of tidal stripping, shocking or evaporation. That extra-tidal material 
(stars or gas clouds) traces the orbit of the satellite or globular cluster has long 
been known to be a powerful probe of the potential of the Galaxy in the halo, 
and has been exploited extensively particularly in the case of the Magellanic 
Clouds and Magellanic Stream (Murai & Fujimoto 1980, Putman et al. 1999, 
this volume). 

Among the future astrometric missions, SIM is a pointed observatory with 
the main science goal of finding nearby planets, while GAIA (under study) will 
continuously monitor the relative positions of about 10 9 stars above ~ 20 mag. 
limit over the full sky over a period of 4 — 5 years with additional information of 
radial velocity to better than 3 — lOkm/s accuracy (depending on spectral type) 
for stars brighter than ~ 16 mag. (cf. the compiled instrument specs in Table 1, 
and the original documents by Gilmore et al. 1998, Unwin et al. 1998). GAIA 
represents an improvement over the Hipparcos mission by a factor of about a 
thousand in precision, more than a million in the probing volume, and a factor 
of ten thousand in number of objects. 

Helmi, Zhao & de Zeeuw (1999, this volume) show that streams can be 
identified by as peaks in the distribution in the angular momentum space, mea- 
surable with GAIA. Once identified, we can fit a stream with an orbit or more 
accurately a simulated stream in a given potential. Johnston, Zhao, Spergel 
& Hernquist (1999) show that a few percent precision in the rotation curve, 
flattening and triaxiality of the halo is reachable by mapping out the proper 
motions (with SIM accuracy) and radial velocities along a tidal stream 20 kpc 
from the Sun. In particular, they show that the fairly large error in distance 
measurements to outer halo stars presents no serious problem since one can 
predict distances theoretically using the known narrow distribution of the an- 
gular momentum or energy along the tails associated with a particular Galactic 
satellite. We expect these results should largely hold for streams detectable by 
GAIA. These numerical simulations are very encouraging since they show that it 
is plausible to a learn great deal about the Galactic potential with even a small 
sample of stream stars from GAIA. Some unaddressed issues include whether 
stream members will still be identifiable in angular momentum in potentials 
without axial symmetry, and the robustness of both methods if the Galactic 
potential evolves in time. 

Here we illustrate how the properties of tidal streams evolve in time-dependent 
potentials and discuss whether members of such streams might still be identified 
using the 6D information from GAIA. We study the effect of contamination from 
field stars, evolution and non-axial symmetry of the potential, and the lack of 
very accurate radial velocity and parallax on our ability to detect streams. We 
concentrate on satellites which fall in and are disrupted recently (about 4 Gyrs 
ago, well after the violent relaxation phase) and maintain a cold spaghetti-like 
structure. We also put the satellites on relatively tight orbits but which lie out- 
side the solar circle (pericenter of about 8 kpc and apocenter of about 40 kpc) 
such that the bright member stars in the stream are still within the reach of 
detectability of GAIA. Such streams typically go around the Galaxy less than 5 
times since disruption, and are typically far from fully phase-mixed. 
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Table 1. Preliminary specifications of SIM and GAIA, as adopted in 
this paper 



program 

launch + duration 
orbit 

resolution (mas) 

parallax (fxas) 

proper motion (fiasyr^ 1 ) 

radial velocity (kms -1 ) 

wavelength(/xm) 

objects (mag. limit) 

coverage 



SIM 

NASA origin 
2006 + 5 

Sun-synchronous (lOOmin) 

~ 10 

4 

1 

0.4-1 

10 4 (V ~ 20) 

pointed, narrow/wide field 



GAIA 

ESA (study) 
2009 + 5 
L2 

~ 100 

4-200 for V < 10 - 20 mag. 
4-200 for V < 10 - 20 mag. 
3-10 down to V ~ 16 mag. 
0.4-1 

10 9 (V ~ 20) 
full sky scan 



2. Strategy for GAIA: stars in a stream vs. field stars 

We propose to select bright horizontal/giant branch (HB/GB) stars as tracers 
of tidal debris of a halo satellite (which we take to be either a dwarf galaxy 
or a globular cluster). A satellite with a typical luminosity L = 1O 5_7 L0 has 
numerous HB and GB stars with My < 0.75mag, which are easily observable at 
20 kpc from the Galactic center (my pa 18mag or brighter) with GAIA. While 
the depth of the debris is typically difficult to resolve with GAIA parallax, GAIA 
proper motion (cf. Table 1) is good enough to resolve the internal dispersion of 
the debris, which is of the order cr/20kpc pa 100/uas/yr, where a ~ lOkrns -1 is 
the typical velocity dispersion of a satellite. 

A satellite with luminosity 10 5 " 7 L Q typically has 

N obs = f L > 100 (1) 
5OOL - v ' 

numbers of stars in its tidal tail, where we assume there is about one HB or 
GB star per 500L Q and between / = 0.5% to 50% of the stars in the original 
satellite are liberated. 

In comparison, the density of field HB or GB stars in the halo which happen 
to be in the same solid angle, and have the same proper motion and radial 
velocity 

N field = Sfi ( r) ~ 1 (2) 

; VlOOkms- 1 / V ' 

where 100 km s -1 is the dispersion of field stars, S ~ (1 — 10) is the typical 
number of field giants per square degree and is the solid angle of the tidal debris 
in degrees, which for the Sagittarius (Sgr) tidal stream is about 5° x (20° — 60°). 
In other words while there might be 25 stars in a solid angle of 5° x 5° in a stream, 
only 25% of the chance that there is one field star in the same piece of sky with 
proper motions and radial velocities indistinguishable from the stream stars. In 
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fact Sgr was discovered on the basis of radial velocity and photometric parallax 
against a dense foreground of bulge stars (Ibata, Gilmore, & Irwin 1994). So as 
far as identifying stars in a cold stream with GAIA is concerned we conclude 
that contamination from field halo stars is likely not a serious problem. 

It is worth commenting on the advantages of stars in a stream, as compared 
to stars in the field, in constraining the Galactic potential. Stars in a stream 
trace a narrow bunch of orbits in the vicinity of that of the center of the mass, 
and are correlated in orbital phase: they can all be traced back to a small volume 
(e.g., near pericenters of the satellite orbit) where they were once bound to the 
satellite. Hence we expect a tight constraint on parameters of the Galactic 
potential and the initial condition of the center of the mass of the satellite 
(about a dozen parameters in total) by fitting the individual proper motions 
of one hundred or more stars along a stream since the fitting problem is over- 
constrained. In contrast, field stars are random samples of the distribution 
function (DF) of the halo, and the large number of degrees of freedom in choosing 
the 6-dimensional DF makes the problem often under-constrained: one generally 
cannot make the assumption that the halo field stars are in steady state as an 
ensemble because it typically takes much longer than a Hubble time to phase-mix 
completely at radii of 30 kpc or more. 

How well can we determine the Galactic potential with a stream? Assume 
each debris star has an intrinsic energy spread AE ~ V? r a/R relative to the 
orbit of the center of the mass, where a = (0.1 — 1) kpc is the size of the parent 
satellite at the time of disruption, also the thickness of a cross-section of the 
stream, and R = (8 — 40) kpc is the radius of the pericenter, then the accuracy 
of the potential from fitting N t, s ~ (100 — 1000) stars is 

1 Af? 

*~-m=W~ (0.1-1)% (3) 

The accuracy is not very sensitive to the luminosity of the disrupted satellite 
since both iV o6s and a scale with the luminosity of the satellite. While this 
agrees fairly well with the study of parameterized static potentials (Johnston et 
al. 1999), it remains to be tested for more realistic models with evolution. 



3. Science for GAIA: evolution of the Galactic potential 

To study the effect of the evolution and flattening of the potential on a stream, 
we follow the disruption of satellites in a simple, flattened, singular isothermal 
potential <3?(r, #,i) which is time-dependent but maintains a rigorously flat ro- 
tation curve at all radii, where (r, 9) are the spherical coordinates describing 
the radius and the angle from the North Galactic Pole, t is defined such that 
t = would be the present epoch. This is perhaps a sensible choice since we 
do not expect a major merger to have occurred within the past 4 Gyrs. We 
study three types of moderately evolving potential with identical flattening and 
rotation curves at the present epoch. 

First we consider a Galactic potential 



A s log r + - cos 29 



(4) 



where 

2irt 

A s (t) = 1 - e + e(t), e(t) = e cos — -. (5) 

Jg 

This model simulates the effect of the Galaxy becoming more massive and flat- 
tened in potential as it grows a disk. The time-evolution is such that the Galactic 
potential grows from prolate at time t = —Tq/2 to spherical at time t = — lb/4, 
and then to oblate at t = 0. A more general prescription of the temporal varia- 
tion might include a full set of Fourier terms. We adopt parameters 

y = 200kms-\ eo = 0.1, (6) 

such that the present-day rotation curve amplitude is 200 km s -1 , and equal 
potential axis ratio ~ 1 — e ~ 0.9. eo needs to be small so as to guarantee 
that the volume density of the model is positive everywhere at all time. We set 
Tg/4 = 4Gyr, a reasonable time scale for the growth of the Galactic disk. 

Second we consider a Galactic potential where the minor axis of the poten- 
tial slowly flips over a time scale Tp, 



logr + — cos 2(6 — 
B 2 v T F ' 



(7) 



This potential might mimic the tidal harassment of Local Group galaxies. We 
set Tp = 2Gyr for the time to flip 180°. This type of model has also been 
examined by Ostriker & Binney (1989). 

For the last case we analyze a static Galactic potential plus a time-varying 
perturbation coming from a massive satellite on a circular orbit with a period 
Tp. The potential is 



<f> P (r,6,t) = V 2 



log r + ^ cos 20 



GMP (8) 



where we use a Plummer model for the perturber's potential with a scale length 
bp about half of the tidal radius, the total mass of the perturber Mp is taken to 
be 5% of the Galactic mass enclosed by its circular orbit of period Tp and radius 
Rp = V °^~ p ■ Generally, the orbital plane of the perturber is unrelated to that 
of the stream, but here they are taken to be the same plane with both moving 
in the same sense. The perturber is set on an exactly polar, circular orbit with 
a period T P = 2Gyr (R P = 64 kpc, M P = 3 x 10 10 M o , 2b P = 14 kpc) without 
considering effects such as dynamical friction and mass-loss. These parameters 
might be appropriate for a very massive perturber such as the progenitor of the 
Magellanic Clouds (Zhao 1998a,b), which could significantly "harass" all smaller 
satellites in the halo. 

Fig. 1 shows the orbit and morphology of the simulated stream in these three 
potentials. The orbit of the disrupted satellite is chosen so that the released 
stream stays in the polar xz plane, which passes through the location of the Sun 
and the Galactic center; the xyz coordinate is defined such that the Sun is at 
x = —8 kpc and y = z = 0. 

We then simulate observations of mock data of 100 bright HB and GB 
stars convolved with GAIA accuracy. The particles in the disrupted satellite are 
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initially distributed with an isotropic Gaussian density and velocity profile (as in 
Helmi Sz White 1999) with dispersions 0.4 kpc and 4kms _1 respectively. These 
particles are released instantaneously at the pericenter 8 kpc from the center 
4 Gyrs ago. These parameters might be most relevant for satellites such as 
the progenitor of the Sagittarius stream. The parallax and radial velocity from 
GAIA will not be very constraining for stream stars beyond 10 kpc (WOfias 
in parallax) because of rapid growth of error bars with magnitude of the stars 
(Lindegren, private communications). Here we adopt a simple parameterization 
for the error of the sample stars. We consider only horizontal branch stars, 
for which the errors in parallax (tt in ftas), proper motion (// in fias yr _1 ) and 
heliocentric radial velocity (Vh in kms -1 ) are functions of the parallax. We find 
that this simple formula 

<7 W = 1.6a,, = oy h = /(tt), f{x) = 5 + 50(50/x) L5 , (9) 

approximates the GAIA specifications fairly well. 

Nevertheless, Fig. 2 and Fig. 3 show that the distribution of the proper mo- 
tions of the stream remains narrow (as shown by the proper motion vs. proper 
motion diagram and the position-proper motion diagram) after taking observa- 
tional error into account. The narrow distribution allows stars in a stream to 
be selected out from random field stars, as we argue in the previous section. 

Fig.4 shows the simulated streams in the energy and angular momentum 
space for various evolution histories of the Galactic potential. For example, in 
the model where the Galaxy slews, particles move along straight lines defined 
by the Jacobian integral E — ^J, where J is the angular momentum in the 
direction around which the system slews. By and large the energy E of particles 
across each stream is spread out only in a narrow range at each epoch in the 
three models; the same holds but to a lesser extent for the angular momentum 
vector J. This implies that stars in the stream are largely coeval even in the 
presence of realistic, moderate evolution of the Galactic potential. The energy 
and angular momentum is also modulated with particle position in a sinusoidal 
way across the stream, an effect which in principle can be used to infer the 
evolution rate of the Galactic potential and the flattening of the potential. 

One of the challenges of using streams to constrain the potential is that 
measurements of both parallax and radial velocity from GAIA are likely domi- 
nated by noise for the fainter [My ~ 20 mag), more distant (~ 50 kpc) members 
of a halo stream (cf. Table 1). Fortunately, we can use the property of nearly 
constant angular momentum and energy across the stream to predict the miss- 
ing information with an accuracy often comparable to or better than directly 
observable by GAIA. In essence we apply variations of the classical method of 
obtaining "secular parallaxes". The simplest example is a polar stream, which 
has no net azimuthal angular momentum, J z ~ 0, so parallaxes can be recov- 
ered (to about lO^as accuracy) from the solar reflex motion of the stream in the 
longitude direction as shown by the linear regression ir ~ |w|/40 in the top right 
panel of Fig. 2. More generally we can use the property that the total angular 
momentum and an approximate energy are roughly constant, i.e., 

J = r x V ~ constant, E ~ (200 kms -1 ) 2 logr + -V 2 ~ constant (10) 
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to predict both parallax and heliocentric velocity; here we simply pretend that 
the Galactic potential is spherical. Surprisingly, these very rough approxima- 
tions often yield fairly accurate parallaxes (~ 10%) and heliocentric velocities 
(~ 30kms _1 ) as shown by the narrow bands in Fig. 2 and Fig. 3; predictions 
tend to be poorer for particles in the (anti-)center direction, because it is where 
the angular momentum J becomes insensitive to the heliocentric velocity. The 
predicted velocities and parallaxes are testable with direct observations at least 
for the nearby brighter members of a stream. 

Most relevant to constraining the evolution of the Galaxy is that the differ- 
ences in the position-proper motion diagram (Fig. 5), reflected by these different 
evolution scenarios, are at a level clearly resolvable with GAIA accuracy. We 
caution that the examples shown here are biased towards nicely-structured and 
extended streams. These happen quite often for models where the Galaxy simply 
slews, but not for models involving a growing disk or a massive perturber as we 
increase the amplitude of the perturbation. We have also run simulations with 
various parameters for the satellite orbit and initial size and the Galactic poten- 
tial. The structure of a stream can become very noisy for highly eccentric orbits 
with a pericenter smaller than 8 kpc and/or for potentials where the temporal 
fluctuation of the rotation curve is greater than 10%. These noisy structures as 
a result of strong evolution can be challenging to detect. Nevertheless we con- 
clude that tidal streams are excellent tracers of the Galactic potential as long 
as a stream maintains a cold spaghetti-like structure, in particular, the results 
of Johnston et al. (1999) and Helmi et al. (1999) for static Galactic potentials 
are likely to be largely generalizable to moderately time-evolving potentials. 
However, perhaps the most exciting implication of these preliminary results is 
that by mapping the proper motions along the debris with SIM or GAIA we 
could eventually set limits on the rate of evolution of the Galactic potential, and 
distinguish among scenarios of Galaxy formation. 



HSZ thanks Amina Helmi and Tim de Zeeuw for many helpful comments 
on an earlier version. 
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Figure 1. Simulated orbits of the parent satellite and the released 
stream of about 100 giant stars in potential models with a growing disc 
&c (upper left panel), a perturber i>p (upper right) and a flipping disc 
(the lower larger panel). The Galactic center is marked by a cross. 
We assume that the parent satellite is totally disrupted 4 Gyrs ago at a 
pericenter 8 kpc (Sun's mirror image point) and moving south (down) 
with a velocity 380 km s -1 . The stubby debris tail shortly (0.3 Gyrs) 
after the disruption is also shown together with the well-developed tidal 




Figure 2. Mock data of the stream convolved with GAIA accuracy 
for the potential model with a flippig disc $_p. The two small panels at 
the top are the proper motion vs. proper motion (jU/ vs. fib) diagram 
and the (longitude) proper motion vs. parallax (/z; vs n) diagram; the 
linear regression \fii\ ~ 200 kms _1 /D kpc ~ AOtt manifests the reflex 
of the Galactic rotation of the Sun and the absence of any azimuthal 
rotation of a polar stream. The lower three panels show the heliocentric 
radial velocity, the proper motion, and the parallax across the stream 
as functions of the latitude; the GAIA error bars for the proper motion 
are typically much smaller than the size of the symbols. The two narrow 
bands bracket the range of parallax^and heliocentric velocity predicted 
from assuming a constant angular momentum |J| and approximate 
energy E across the stream; the width of the bands here are computed 
from the initial spread of E and J. 
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Figure 4. The lower panels compare the angular momentum and en- 
ergy along three streams simulated in three different time-varying po- 
tentials. Each shaded region in the top panel shows how most (~ 70%) 
stars in a stream evolve in the energy vs. angular momentum plane for 
one of the three potentials. They are shown in steps of about 0.15 Gyr 
from 4 Gyrs ago to the present epoch (marked by the symbols). The 
plus symbols are for the flipping Galactic potential with Tp = 2 Gyrs, 
the diamonds are for the case with a massive perturber on a Tp = 2 
Gyrs orbit, the crosses are for a slowly growing potential with Tq/A = 4 
Gyrs. J c and E c are the present values for the center of mass of the 
stream. 
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Figure 5. Mock data for the norm of the proper motion vector across 
three streams simulated with GAIA accuracy. The three streams are 
generated in three different evolution histories of the Galactic potential. 
Symbols have the same meaning as in the previous figure. 
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